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I apply a recently developed algorithm for reweighting simulations of lattice QCD from one quark 
mass to another to simulations performed with overlap fermions in the epsilon regime. I test it 
by computing the condensate from distributions of the low lying eigenvalues of the Dirac operator. 
Results seem favorable. 



I. RECALLING THE ALGORITHM 

OO ■ Several recent papers [l|, 0] have described the idea of performing simulations of QCD at one value of the quark 
. mass by reweighting the configurations generated in a simulation at another quark mass. Overlap fermionsjj, 0], 
^vq ■ the most easy to analyze and and most expensive to simulate of all lattice discretizations of the Dirac operator, are 
natural candidates for reweighting. I describe some rather trivial numerical experiments I have made, reweighting 
simulations with two flavors of overlap fermions. The techniques are basically those of A. Hasenfratz, R. Hoffmann 



and S. Schaefer but the properties of the overlap operator eliminate many technical difficulties. 
. Let's summarize the implementation of reweighting I will follow, by considering a single flavor of fermion with Dirac 
' operator D{m). We will not worry about questions of positivity or reality for the moment. Let's suppose that we 
have generated a stream of configurations with one value of the quark mass mi and want, by reweighting, to use the 
stream to perform simulations at another quark mass m2- To do so, one must rcweight configurations by a factor 
p i' 

D: deti^(m2) 
, ' Introducing the operator 

\0- n = D{m2)-'D{mi) (2) 

• we can write w as an average over a set of complex random vectors ^, 

o ■ 

00 . 

' With the usual definition of the link gauge variable U and pure gauge action Sq, expectation values of operators 
J> , evaluated at mass m2 are given by 

X 
J3 



{0)ra, = ^ J dUe-^^detD{m2) 0{U) 
Z2 J 

= |i(O(t/)exp(-e^(f}-l)0>™„5 
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(4) 



where 



^ = (cxp(-Ct(f7-l)0)™,.C. (5) 

If we imagine that our simulation at mass mi consists of a stream of pairs of variables {!7i,^i}, then the expectation 
value is 

E.exp(-eJ(J7(t/0-l)?0 
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That is, we reweight each configuration by a factor 

= eM-CKmd - m)- (7) 

This was so far aU quite generaL Now we assume that we are doing simulations with overlap fcrmions. For any 
number of flavors, all calculations can be performed using the squared Hermitian Dirac operator D{miy D{mi) = 
H^mf) — SiH{0)'^ + mf where Si = 1 — mf /(4i?o)^ and H{0) is the massless squared overlap Hermitian Dirac operator 

i/(0)=ro(75-eW). (8) 

The quantity h is the kernel operator h — ^^{d — Rq) (in terms of a kernel Dirac operator d) and e{h) is the matrix 
sign function. We assume that we have recorded a set of eigenfunctions of H{0) (and their associated eigenvalues), 
H{0)\k) = Xk\k). The spectrum of iJ(0)^ consists of a set of zero eigenvalue chiral modes and a set of degenerate 
(paired) nonzero eigenvalue eigenmodes of opposite parity. The nonzero mode contribution to w in Eq. [3] can be 
computed using random vectors ^ which are chiral, with chirality in the sector without zero modes. Each flavor of 
dynamical fermion has its own chiral random vector. 

Now we come to the question of practicality: Reweighting will fail if the weight of each configuration deviates 
widely from the mean, because then only the (presumably small number) of configurations carrying a large weight 
will contribute to averages. It can also fail if the estimator (Eq. [7|) has a large variance, for then one will need to 
average the same underlying gauge configuration over many estimators. Can schemes be devised, so that the weights 
Wi do not fluctuate too much from configuration to configuration? Presumably what will work will depend on the 
simulation and reweighted quark mass and the simulation volume. The phase space of possible choices is large. 

It is always a good thing to replace as much of the stochastic estimator of the determinant with an exact result. 
Introducing the Hermitian projector onto low eigenmodes of _ff(0)^ (call it P and its complement _P = 1 — P), 

detn = detPn detPn. (9) 

We compute Pf2 from eigenvalues and we only need to make a stochastic estimator for the high eigenmode part of 
the weight Wi 

^ \2 I 2 

w. = i^f'' -Q '^^k + ^2 exp(-et(TOP - m ^ wi^Mgh (10) 

(reweighting a configuration with A^o zero modes from mass mi to m2, and considering a single fiavor). To complete 
the equation set, 

where ci2 = m\ — (si/s2)"i2, and with y = because the random vector can only live in the space of PftP, 

Whigh = expi-iy^ri- l)y). (12) 

Projection of low modes plus the use of a random vector in the chirality space without zero modes improves the 
effective conditioning number of H(m2)^^ ■ 

One interesting place in parameter space is the epsilon regime. Here the quarks are so light that the pion "fills 
the box" - if the volume is = L'^, then mjrP << 1 (and all other mass scales M large, ML > 1) defines the 
epsilon regime. Let's perform some experiments there: I have several sample data sets with Nf = 2 flavors of overlap 
fermions on 12^ simulation volumes at a nominal lattice spacing of 0.14 fm. I will use a quark mass aniq ~ 0.03 
fm (nominally about 43 MeV) and am, = 0.01. They were generated using the hybrid Monte Carlo algorithm, with 
the reflection/refraction algorithm devised in Ref. [5]. They used the differentiable hypercubic smeared link of Ref. 
Q and one or two additional heavy pseudo-fermion fields as suggested by Hasenbuschjv | . The integration is done 
with multiple-time scalesjsj. Details of the actions are described in Refs. 10, 11, [l2, [l^. All in all, these are 



very conventional overlap fermions. I typically compute the lowest 12 eigenvectors and eigenvalues of H{0)^; these 
eigenvalues run up to about A ~ 0.04. 
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FIG. 1: Determinant weight from a set of aniq = 0.03 configurations into (a) aiUq = 0.01, (b) am, = 0.035, normalized by the 
average weight, computed without removing eigenmodes. Panels (c) and (d) show the reweighting to mass aniq = 0.01 and 
0.035 while treating the 12 lowest eigenvalues of H{0) exactly. 



The first test checks what can gained by removing eigenmodes. We take a set of 22 lattices from our stream of 
amq = 0.03 simulations and reweight them to a set of target masses: aniq = 0.01 and 0.035. In this test we averaged 
the stochastic part of the weight over six pairs of two chiral pseudofermions. Fig. [T]show a comparison of the resulting 
weights either not keeping any eigenmodes or removing the lowest 12 eigenmodes from the stochastic estimator. The 
error bars show the variation in weight over the ensemble of pseudofermion noise vectors used for each configuration. 
It is clear that the choice of removing eigenmodes is the superior one, from the point of view of suppressing the 
variance of the estimator. 

The low eigenvalues do not capture the entire reweighting factor. Fig. [2] shows the weights from just the low 
eigenvalues (divided by the average reweighting factor from the true weights). Their (incorrectly normalized) values 
appear to track to full weight. 



II. EIGENMODE DISTRIBUTIONS IN THE EPSILON REGIME 

By themselves, pictures of the fluctuating weights give no indication of how well an actual reweighted calculation 
will perform. A test is needed. For a little physics example I select the problem of determining the condensate 
from the distribution of low- lying eigenvalues of the massless Dirac operator in the epsilon regime in sectors of fixed 
topology. These distributions are given by Random Matrix Theory (RMT) 14 , 1^, [3, Hi ■ Overlap fermions are 



optimal for this project (as for any epsilon regime simulation) due to the control they give over lattice topology. 

I am aware of three previous measurements of from eigenvalues with Nf = 2 flavors of dynamical overlap 
fermions. Two of them, Refs. [12 and [13. were not really in the epsilon regime; in the second paper, the bare 
quark mass is amq = 0.03 corresponding to a pion mass in lattice units of arriTr = 0.324 (so rriTri ~ 3.9). The JLQCD 
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FIG. 2: Contribution to determinant weight from the 12 lowest eigenvalues of H{0)^ for the am, — 0.03 configurations of Fig. 
□ (a) amg = 0.01, (b) am, = 0.035. 



collaboration, Ref.[19], has a true epsilon regime calculation of E. I am also aware of two recent studies which use 
dynamical fermions which are not exactly chiral, but which are said to have highly improved chiral symmetry: Refs 
[20| and 21|. (The latter simulation used 2+1 flavors.) All of these papers produce similar and unsurprising values 
for the condensate, S ~ (250 MeV)"^. An unpleasant feature of the epsilon regime is that finite volume corrections 
are power law, not exponential. The effect is to replace the value of the condensate extracted from the RMT fit, E, 
by T,L — psE where 



iV? - 1 1 



(13) 



with A(0) the contribution to the tadpole graph (propagator at zero separation) from finite-volume image terms. In 
the epsilon regime, A(0) = —f3o/VV and /3o depends on the geometry 2 2|. (It is 0.1405 for hypercubes.) 

I carry two data sets into - or closer to - the epsilon regime. The first data set is the set of arriq — 0.03 configurations 
from Ref. a stream of 400 thermalized hybrid Monte Carlo trajectories of unit length giving 30 = lattices 

and 75 \v\ = 1 ones. The second set is a collection oi 7A ly — lattices from four little streams (each of about 150 
trajectories) run at amq — 0.01. This is already on the edge of the epsilon regime. Nearly all were f = 0. I reweighted 
the aruq = 0.03 data set to target masses amq = 0.01, 0.015, 0.02, 0.025, and 0.035. I reweighted the aniq = 0.01 set 
to 0.0025, 0.005, 0.0075, 0.025 and 0.03. I used only a single pair of pseudofermions per configuration, since I treated 
the lowest 12 eigenvalues of H{0)^ exactly. 

The data were analyzed with a conventional bootstrap analysis. In the bootstrap, the weight of a configuration 
was the number of times it was selected for the bootstrap, times the (normalized) weight factor from the determinant 
ratio. Results for Y^lV from a fit to the lowest eigenvalue distribution in one topological sector are shown in Fig. [3l 
Of course, all the reweighted points at different quark masses are highly correlated; they came from the same data 
sets. It appears that reweighting into the epsilon regime was successful, while trying to go to larger quark masses 
(arUq ~ 0.01 to 0.035, for example) was less so. It is probably no surprise that reweighting for a small change in mass 
works better than reweighting a big change. 

Readers might recall that to complete a calculation of El, one needs a separate determination of a lattice spacing 
and a lattice-to-continuum matching factor Zs- Zs was determined for this action in Ref. [3]: Zg^^{2 GeV)=0.76(3). 
Of course, the lattice spacing varies as the bare parameters of the simulation change. However, this variation is small 
in the epsilon regime simply because the absolute change in the quark mass is small. For this data set, ro/a = 3.71(5) 
at aniq = 0.03 and 3.77(7) at aniq = 0.01. Thus the three unreweighted values of TqE^ from this study are 0.326(30) 
{amq = 0.01, ly = 0), 0.347(37) {amq = 0.03, = 0), and 0.294(20) (am, = 0.03, = 1, and with ro - 0.5 fm. 
El - (260- 270 MeV)^. 
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FIG. 3: The quantity I^lV from reweighted 12* lattices, as a function of the target reweighted mass, (a) True mass aruq = 0.01, 



u = 0. Crosses are reweighted points, the octagon is unweighted, (b) True mass aniq = 0.03: u 
V = 1 lattices are marked with squares when reweighted, diamond when not. 



lattices labeled as in (a) 



III. CONCLUSIONS 



Reweighting dynamical overlap fermion data sets into the epsilon regime worked better than I expected. Groups 
doing simulations with overlap fermions might well be advised to investigate it as a technique. All the ingredients will 
probably already be in hand. The main reason for reweighting nonchiral actions - namely, that one wants to avoid 
exceptional configurations - obviously does not apply to overlap fermions. However, overlap fermion simulations are 
so expensive that running at many parameter values is daunting. Any methodology which allows one to recycle old 
configurations is worth exploiting. 
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